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We use the multifractal formalism to describe the effects of dissipation on Lagrangian velocity 
statistics in turbulent flows. We analyze high Reynolds number experiments and direct numerical 
simulation (DNS) data. We show that this approach reproduces the shape evolution of velocity 
increment probability density functions (PDF) from Gaussian to stretched exponentials as the time 
lag decreases from integral to dissipative time scales. A quantitative understanding of the departure 
from scaling exhibited by the magnitude cumulants, early in the inertial range, is obtained with a 
free parameter function D(h) which plays the role of the singularity spectrum in the asymptotic 
limit of infinite Reynolds number. We observe that numerical and experimental data are accurately 
described by a unique quadratic D(h) spectrum which is found to extend from h m in ~ 0.18 to 
hmax ~ 1, as the signature of the highly intermittent nature of Lagrangian velocity fluctuations. 
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Statistical properties of homogeneous three dimen- 
sional turbulence have been studied for a long time in 
the Eulcrian framework Q. Recently a growing inter- 
est in studying intermittency from a dynamical point of 
view has been motivated by high precision Lagrangian 
experiments. Essentially two experimental groups have 
performed particle tracking in highly turbulent flows. 
The group at Cornell |2| reports measurements of La- 
grangian acceleration in a turbulent water flow between 
two counter-rotating disks for Taylor-based Reynolds 
numbers 200 < R\ < 900. The experiment carried 
out at ENS-Lyon in a similar von Karman flow, is 
based on acoustic tracking. It provides Lagrangian veloc- 
ity records covering the inertial range of turbulent mo- 
tion, up to several integral time scales. In addition to 
these complementary experiments, DNS of the Navier- 
Stokes equations 0, 3 have produced comparative nu- 
merical results in the range 75 < R\ < 380. Both ex- 
perimental and numerical studies have revealed the ex- 
istence of a very strong intermittency in the Lagrangian 
dynamics: the PDF of the velocity is Gaussian, while the 
PDF of the acceleration exhibits extremely large tails. 
In between, the PDFs of velocity increments changes 
continuoulsy between these two functional forms as the 
time lag is decreased from integral times (the Lagrangian 
auto-correlation time) to dissipative ones (below the Kol- 
mogorov time scale). Several stochastic models have been 
proposed that reproduce the behavior of the Lagrangian 
acceleration p| and of the Lagrangian velocity increments 
@. They rely on different physical assumptions and all 
include several ad-hoc hypothesis and parameters in order 
to fit the experimental observations — a rather compre- 
hensive review can be found in Q • The aim of the present 
work is not to add to these models but to provide a com- 
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prehensive description of the Lagrangian intermittency, 
using a formalism that describes both the inertial and dis- 
sipative range of time scales. It is motivated by the desire 
to analyze globally the Lagrangian scales, which are likely 
to be dynamically connected (for example it is shown in 
[j| that the Lagrangian acceleration is influenced by the 
large scale dynamics, a feature that is also pointed out 
by models based on rapid distortion ideas |9(). In ad- 
dition we require that the description of intermittency 
be independent on model assumptions. The multifractal 
description jlfj, already widely used in Eulerian studies 
of turbulence, is a natural choice. To encompass iner- 
tial and dissipative features, we recast the multifractal 
picture of the intermediate dissipative ran ge, originally 
proposed for Eulerian velocity fluctuations [Tl], [T^] , into 
the context of Lagrangian velocity. The resulting anal- 
ysis provides a synthetic and comprehensive description 
of the various sets of experimental and numerical data. 

In the present description, a first-order Lagrangian ve- 
locity increment over a time scale r is written as : 

S T v(t) = v(t + t) - v(t) = (3 (r/T) S T v , (1) 

where all the time scale dependence is contained in the 
independent random function (3{t/T) (> since La- 
grangian velocity increment PDFs are symmetric), T be- 
ing the integral (Lagrangian) time scale above which ve- 
locity increments become uncorrelated. The PDF of in- 
tegral time scale increments Stv is thus assumed to be 
Gaussian (Q) — a result of a central limit argument, also 
in agreement with Eulerian observations. Once the dis- 
tribution of V((3) is known, the PDF of increments at 
any time scale r is computed as 

n^-ffofflw. (2 , 

In the standard mutifractal formalism |10| , (3 is assumed 
to have a power law scale dependence in the inertial 
range, f3 ~ (r/T)' 1 , with a spectrum V(h) (meaning that 
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the PDF of observing an exponent h at scale r is pro- 
portional to (t/T) 1 ~' d ^). In the inviscid limit, h and 
T>{h) acquire the mathematical status of Holder exponent 
and singularity spectrum, respectively. Note that this 
description based on first order increments is restricted 
to exponents h < 1 |l3j |. a limitation which will be ad- 
dressed later. To describe the entire range of scales cov- 
ered in experimental measurements and computer sim- 
ulations, one must take into account the effects of vis- 
cosity (finite R\). In the dissipative range, velocity fluc- 
tuations are smoothed by viscous damping (or by mea- 
surement filtering) and the velocity increments become 
proportional to the time scale 5 T v(t) = ra(t), where a(t) 
is the Lagrangian acceleration. Wc shall consider that 
the cross-over between inertial and dissipative statistics 
occurs when the local Reynolds number is of order unity, 
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where Re is the integral scale Reynolds number. This 
defines a local Kolmogorov dissipative time T v (h) = 

TRe 2 ^ 1 where the local velocity increments change from 
inertial scale invariance to dissipative scaling (this im- 
plies h > —1/2, for time scales to be shorter than T). 
Changing the integration variable from /3 to h in Eq. (2) , 
the PDF of velocity increments at any scale t/T can be 
written as the sum of two contributions 

P(M = L /2 dh fk&h) 9 [-Km J 
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where the functions ftd and Vi^d have the proper in- 
ertial (ft - (r/T) h , Vi ~ (t/T) 1 -^) and dissipa- 
tive (ft ~ t/T, V d ~ {T n {h)/T)^ v ^) scalings. The 
change occurs at the critical value h* for which the local 
Reynolds number is unity: 
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For h < h* (t /T, R e ), the increments 8 T v are in the in- 
ertial range, while they lie in the dissipative range for 
h > h* . Finally, we impose that the function j3(r/T) 
be continuous and differentiable at the transition, fol- 
lowing a strategy used in the Eulerian domain [12|, and 
inspired from an elegant interpolation formula proposed 
by Batchelor 0. In this framework, a single function 
(3(h; t/T, Re) covers the entire range of scale 
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with probability density function 
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The PDFs of Lagrangian velocity increments are then 
computed using Eqs.(4-7), and compared to the experi- 
mental / numerical data. The entire range of time scales 
is covered when the four parameters T, Re, 8 and T>(h) 
are prescribed. T and Re, both global parameters im- 
posed by the flow production, are explicitly incorporated. 
T>(h) is a parameter function which must be derived from 
the experimental or numerical data. In this work, and 
as a posteriori justified, we assume a quadratic form 
T>(h) = 1 — (h — ci) 2 /2c2 (a lower order approximation 
customary in intermittency studies). The constants C\ 
and C2 are constrained to satisfy c\ = 1/2 + C2 in or- 
der for Kolmogorov scaling to be observed in the inertial 
range ((S T v 2 ) oc r). S is a free parameter that accounts 
for the smoothing of the transition from inertial to dissi- 
pative scales. In his mode ling of Eulerian velocity struc- 
ture function, Batchelor Qj| sets 8 = 2, a value used 
in other Eulerian studies [12j and supported by theoret- 
ical arguments 15]. In practice T and Re are computed 
using the experimental data (up to some empirical mul- 
tiplicative constants), and the free parameters C2 and 8 
are estimated using a least square minimization scheme. 

In FigQIa) we show the results for the experimental 
data from the ENS-Lyon group (R x = 740) Q. The ve- 
locity increment PDFs are well reproduced at all scales 
for the set of parameter values 8 = 1.08, c\ — 0.575 and 
C2 = 0.075. This agreement is emphasized in Fig^c), 
where the fourth order moment (8 T v) 4: P(S T v) is shown 
to peak at the same values as the experimental curves 
and this for all considered time scales. In Figsma,c) 
we have also included the Cornell acceleration data at 
Rx = 690 2]. What is remarkable is that with quite 
consistent parameter values, namely 8 — 1.3 and still a 
parabolic T>(h) curve with c\ — 0.579 and C2 = 0.079, 
one again reproduces the experimental acceleration PDF 
with a great accuracy in the tails; in particular, the flat- 
ness F = (a 4 )/(a 2 ) 2 = 56.1 is very close to the exper- 
imental value F = 55 ± 8 2]. For our DNS data at 
Rx = 140 (FigsQJbjd)), the agreement is even more pro- 
nounced than for the experimental data. In particular, 
it is hard to distinguish the fit from the data points in 
Fig^d) . Note that the parameter values obtained to get 
the best fit of the numerical data are not significantly dif- 
ferent from those used in Figsnja,c): T>(h) corresponds 
to ci = 0.586 and C2 = 0.086, while 8 = 1.98 is now quite 
close to the Batchelor value 8 = 2 ^j. The fact that 
(5 is smaller in experimental measurements is due to the 
filtering induced by the finite size of the tracer particle 
and signal processing algorithms. 

A first finding of our analysis is that almost identi- 
cal functions T>{h) are obtained for the three sets of 
data (Cornell, Lyon and DNS), although they cover a 
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FIG. 1: Comparison of the experimental (a,c) and numerical 
(b,d) data for the normalized velocity increment PDF V(S T v), 
where S T v = S T v/((8 T v) 2 ) 5 , with the predictions of the mul- 
tifractal description. («)(a,c): ENS-Lyon experimental data, 
for time lags t/T = 1, 0.35, 0.16 and 0.07, from bottom to top; 
the solid lines are the model fit with C2 = 0.075 and 8 — 1.08. 
(•)(b,d) DNS data calculated for t/T = 1, 0.25, 0.17, 0.11 
and 0.05, from bottom to top; the solid lines correspond to 
parameter values C2 = 0.086 and 8 = 1.98. (o)(a,c) Cornell 
acceleration data, the solid lines are the model predictions for 
C2 = 0.079 and 8 — 1.3. The curves are displayed with an 
arbitrary vertical shift for clarity, and the original <5«-axis for 
the acceleration PDF (o) has been shrunk by a factor of 4. 



wide range of scales and of turbulent Reynolds num- 
bers: the symbols in Fig|2]are undistiguishable, certainly 
within error bars. This suggests that such a parabolic 
T>(h) curve should also be relevant in the limit of in- 
finite Reynolds number. A parabolic singularity spec- 
trum is the hallmark of log- normal statistics as originally 
proposed by Kolmogorov and Obukhov [lfj for Eulerian 
velocity statistics. Our findings for Lagrangian veloc- 
ity statistics is a parabolic T>(h) spectrum centered at 
Ci = 0.58 ± 0.01, significantly larger than the K41 value 
1/2, and of width C2 = 0.08 ± 0.01 (commonly called the 
intermittency exponent). This is significantly larger than 
the corresponding value cf = 0.025 ± 0.003 derived for 
Eulerian velocity data • ft corroborates the fact that 
Lagrangian velocity statistics are more intermittent than 
Eulerian velocity statistics 0. Let us remark that the 
ratio of the Lagrangian to Eulerian intermittency expo- 
nents c 2 /cf = 0.08/0.025 = 3.2 ± 0.2 is very close to 
the value (3/2) 3 = 3.375 which can be predicted using a 
Kolmogorov-Richardson argument 0, 0] . 

To emphasize further the quality of the proposed de- 




FIG. 2: V(h) curves extracted from: (□) ENS-Lyon velocity 
data, (V) Cornell acceleration data and (o) DNS numerical 
data. Are also represented for comparison the Eulerian log- 
normal (thick solide line) and log-Poisson (thick dashed line) 
spectra as well as their Lagrangian counter-parts computed 
using Eq.(9) (solid and dashed lines respectively). 



scription of Lagrangian intermittency, we now consider 
the evolution across the scales of the moments of the ve- 
locity increments, the so-called structure functions 
S(p,r) =< \S T v\ p >. As advocated in Ref.[I3 for Eule- 
rian velocity data analysis, the magnitude cumulant anal- 
ysis provides a more reliable alternative to the structure 
function method. Indeed, it is straightforward to derive 
the relationship between the moments of 7 , (|<5 T w|) and 
the cumulants C„(t) of "P(ln |<5 T v|): 



< \6 T v\" >=exp ^C»(r)^ 



(8) 



In the inertial range, where multifractal power law scal- 
ing is expected to be observed for the structure functions, 
the cumulants should behave like C„(r) lnr. In 

FigO we report the results of the computation of the cu- 
mulant C„(t) for both the ENS-Lyon experimental data 
(FigsEHa-c)) and our DNS numerical data (FigsOHd-f)). 
We have explicitly subtracted the cumulant of the Stv 
Gaussian PDF and the curves are computed using the 
same parameter set as in Fig^ For both the experi- 
mental and numerical cumulants, it is quite convincing 
that our multifractal decription provides a comprehen- 
sive understanding of the observed departure from scal- 
ing when going from large t/T ~ 1 to small t/T ~ 10~ 2 
time scales. In fact, an unambiguous inertial range is ob- 
served for the experimental second order cumulant only 
(FigEJJb)). For the first order cumulant, both the exper- 
imental (Fig|3Ja)) and numerical (Fig[3{d)) curves dis- 
play some curvature over the entire range of scale, a fea- 
ture that is well reproduced by our extended multifrac- 
tal description. For C-zir/T), the model describes quite 
well the cross-over behavior observed in the experimen- 
tal (Fig|3fb)) and numerical (Fig|3{e)) data down to time 
scales t/T of the order of the smoothing filtering scale 
(including the finite size of the tracer particles). The pre- 
dicted plateau C 2 (r/T) -> C 2 (R\) in the limit t/T -> 0, 
is reached in the numerical but not in the experimental 
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sipative range should now scale as (r/T) 2 and (ii) h* may 
be larger than one and the integrals must be computed 
in the order of increasing h's. When using the same pa- 
rameters as before, we observe in Fig|3] that the exper- 
imental and numerical cumulants are robustly and even 
better reproduced. Let us point out that as seen in Figs. 
Etb,c,e,f), because the C„(r/T) for n > 2 are predicted 
to reach a plateau C n (Rx,N) - (1 + 2N) n \nR x in the 
limit of vanishing r's, then the experimental and numer- 
ical cumulants for the second-order increments are much 
more affected by the smoothing filtering process in the 
intermediate dissipative range than previously observed 
for the first-order increments. 

To conclude, we return to our observation that a 
unique T>{h) spectrum yields an accurate description of 
the Lagrangian velocity statistics at all scales. Such a 
spectrum V E (h) has been extensively studied in the Eu- 
lerian domain Two widely used forms (correspond- 
ing to log-normal and log-Poisson statistics) are shown in 
Fig. 2. They can be mapped into the Lagrangian domain 



FIG. 3: Cumulant analysis C n (r/T) vs In r/T for: (a-c) 
ENS-Lyon experimental data and (d-f) DNS data, (o) first- 
order velocity increments (N = 1); (□) second-order velocity 
increments (N = 2). The associated multifractal descriptions 
for N — 1 (solid line) and N = 2 (dashed line) correspond to 
the same parameter values as in Fig.l. Curves in (a,d) are 
vertically shifted for clarity. 

data. For C^t/T), we note that the DNS data in the 
inertial range are quite well reproduced taking C3 = 0, at 
odds with the experiment. 

We now return to the assumption that the scaling ex- 
ponent h should remain smaller than one. It may be seen 
in Fig. 2 that the largest measured h values reach 1. In 
this case, it is known that first order increments are not 
adapted to detect singularities with exponents h> 1 ^| • 

We have thus repeated our analysis using second order 
(2) 

velocity increments S T v(t) = v(t+r) — 2v(t+T/2)+v(t). 
Eqs. (4,6,7) must be slightly modified because (i) the dis- 



V(h) = -h+(l + h)V E (h/(l + h)) , (9) 

using a Kolmogorv Refined Similarity argument in the 
spirit of the work done by Borgas [18| . The resulting 
curves are shown in Fig. 2; we note that the agreement 
with the measured Lagrangian T>(h) functions is excellent 
on the left-hand side of the curves, i.e. for values h < c\ 
corresponding to intense velocity increments. On the 
right-hand side (h > Ci) there is a noticeable difference. 
Whether this difference is significant deserves more inves- 
tigation. It may be of importance since the above rela- 
tionship clearly shows that the Eulerien and Lagrangian 
singularity spectra cannot be both log-normal. 

Numerical simulations were performed at CINES 
(France) using an IBM SP computer. We wish to ac- 
knowledge Bernard Castaing for his critical comments. 
We are grateful to the Cornell team for making their 
data available. 
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